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Abstract 

The B dccaj have been observed for the first time. The branch- 

ing fractions of the B K~ K^*'>^ decay modes are measured. Significant signals 

are found for the B D^*^K~K*^ and B~ D^K~Kg decay modes. The invari- 
ant mass and polarization distributions for the K~K*^ and K~ subsystems have 
been studied. For the K~K*^ sybsystem these distributions agree well with those 
expected for two-body B D^*'^a^ (1260) decays, with a5"(1260) R- K*^ . The 
analysis was done using 29.4 fb~^ of data collected with the Belle detector at the 
e"^e~ asymmetric collider KEKB. 

Key words: B decay, branching fraction, ss pair 
PACS: 13.25.Hw, 14.40.Cs, 14.40.Nd 
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1 Introduction 



In this paper we report on tlie first observation of B meson decays to 
Di*)K-K^*^^ final state^. Sucli decays require tlie creation of an additional 
ss pair and can occur via tfie quasi-two-body meclianism sliown in Fig. la, 
wliere tfie intermediate resonance decays to K~K^^*\ or by a non-resonant 
tliree-body decay sliown in Fig. lb. Contributions from color-suppressed dia- 
grams are expected to be small [1]. The investigation of these processes can 
provide important information for testing various models of hadronic B decays 
as well as for studies of K~K^^*^ resonant states. 




Fig. 1. The external spectator diagrams for a) quasi-two-body and b) non-resonant 
three-body B meson decays with the ss pair creation. The decay channels 
^0 ^ D+K'K*^ and B' D^R-R^ are used as examples. 

In B ^ D^*^ K^*^^ decays, any intermediate K~K^^*^ resonance must have 
isospin 1. The allowed quantum numbers are J^= 0~,1~,1"'' for the K~K*^ 
final state and 1~ for the K~Kg final state. A K~K^ state with J^= 0^ 
cannot be produced via tree diagrams of Figs. la,b [2] in the limit of exact 
isospin symmetry. The production of resonances with spin larger than one is 
forbidden in the factorization approach [1]. Furthermore, these restrictions on 
the quantum numbers for K~K'^^*^ resonant states also hold for non-resonant 
K-K°(*^ production (Fig. lb) [2]. 

In the case of decays via intermediate resonances, the final state branch- 
ing fraction depends on both the production and decay rates of the reso- 
nance. Large production rates of the order of 1% have been observed for the 
B L)(*)ai(1260) and B D^*^p{770) decay channels [3]. The R-R*^ final 
state can be produced in the decay of the aj"(1260) resonance, however, the 
branching fraction for aj"(1260) R~R*^, which was indirectly determined in 
r decays by the CLEO collaboration [4], was found to be only (3.3 ± 0.5)%. 
One also expects that for the J^=l~ resonances p(1450) and p(1700) the 
branching fractions to R~R^^*'> are at the level of a few percent or less [5]. 

^ Charge conjugate modes are implicitly included. 
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For K production in p(770) meson decay the phase space is very hmited. 

In contrast to two-body decays, non-rcsonant three-body B decays to one 
charmed and two hght mesons have not yet been observed [3]. It should be 
noted that ss pair creation in a non-resonant three-body B decay (Fig. lb) 
is only moderately suppressed relative to that for uu and dd pairs if SU(3) 
and phase space effects are taken into accoimt [6]. Unfortunately, quasi-two- 
body and three-body processes cannot be easily separated because of the large 
widths of the intermediate resonances. 

The branching fractions for D decays with ss creation are measured to be 
~ (1 — 2)% [5]. However, these values cannot be directly extrapolated to B 
decay because final state interactions, which are much larger for D decays, can 
significantly affect the decay rates. A possible similarity between the diagrams 
of Fig. 1 and those describing r lepton decays with K'~K^^*'> in the final state 
is obscured by the rather limited phase space, resulting in small branching 
fractions for r decays to K~K^v-,- and K~K*^i/t- when compared to non-ss 
channels. 

In this paper the branching fractions for B — > decays are ob- 

tained for the full allowed kinematic region. The invariant mass and polariza- 
tion distributions of the K~K'^^*^ subsystem are then studied in detail. 



2 Belle Detector 



The data were collected with the Belle detector at KEKB, an asymmetric 
energy double storage ring colhder with 8 GeV electrons and 3.5 GeV positrons 
[7]. The results are based upon a data sample with an integrated luminosity 
of 29.4 fb~^ taken at the T(4S) resonance which corresponds to 31.9 million 
BB pairs. 

Belle is a general-purpose detector with a 1.5 T superconducting solenoid mag- 
net. Charged particle tracking is provided by a silicon vertex detector (SVD), 
consisting of three nearly cylindrical layers of double-sided silicon strip de- 
tectors, and a 50-layer central drift chamber (CDC). Particle identification 
is accomplished by combining the information from silica aerogel Cerenkov 
counters (ACC) and a time-of-fiight counter system (TOF) with specific ion- 
ization (dE/dx) measurements in the CDC. The combined response of the 
three systems provides at least 2.5(J equivalent K/n separation for laboratory 
momenta up to 3.5GeV/c. A CsI(Tl) electromagnetic calorimeter (ECL) lo- 
cated inside the solenoid coil is used for detection of photons and electrons. 
The detector is described in detail elsewhere [8] . 
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The Monte Carlo (MC) simulation event samples used in this analysis were 
generated using a detailed GEANT-based simulation of the Belle detector 
response. Simulated events were processed in a manner similar to the data. 



3 Selection Criteria 

The charged track momenta are reconstructed using the CDC and SVD track- 
ing system. Kaon and pion mass hypotheses are assigned to charged tracks 
using a likelihood ratio P{K/tt) = L{K) / {L{K) + L{n)) , which ranges between 
and 1. P{K/7r) > 0.6 is required for kaon candidates and P{K/t^) < 0.6 for 
pion candidates. 

Candidate 7r° mesons are reconstructed from pairs of photons in the ECL 
with invariant masses within ±15MeV/c^ of the nominal 7r° mass (~ 3(j in 
the 7r° mass resolution) . The 7r° daughter photons are required to have energies 
greater than 30 MeV. A mass-constrained kinematic fit is performed on the 7r° 
candidates to improve their energy resolution. 

Kg candidates are formed from tt'^tt" combinations with an invariant mass 
within ±10MeV/c^ of the nominal mass (~ 3(t). The two pions are re- 
quired to have a common vertex that is displaced from the interaction point by 
at least 0.4 cm in the plane perpendicular to the beam direction. The difference 
in z coordinates (parallel to the beam direction) for the tracks constituting 
the secondary vertex is required to be less than 2 cm. The angle a between 
the flight direction and the measured Kg momentum direction is required 
to satisfy cos a > 0.8. 

An opposite sign K and tt meson combination is taken as a K*^ candidate if 
its invariant mass lies within a ±50 MeV/c^ interval of the nominal K*'^ mass. 

The decay channels — > K~tt^ti^ and K^n^ are used in this analy- 

sis for all studied B decay modes, in order to avoid large combinatorial back- 
ground. The decay channel D° — > ir~7r~7r+7r+ is only added for the D*'^ and 
meson reconstruction, because the D* constraint strongly suppresses the 
combinatorial background. The invariant mass of the D candidates is required 
to lie within ±15MeV/c^ of the nominal D mass for the first two modes and 
within ±12 MeV/c^ for the K~7i~n~^7T~^ mode (~ 3a). A mass and vertex con- 
strained kinematic fit is then performed on the D candidates and results of 
good quality fits are used to improve the D momentum resolution. 

For and D*+ candidates, the D*^ D^tt^ and D*+ D°7r+ decay 
modes are used. The invariant mass of the D'^ and tt combination is required 
to be within ±3.0MeV/c^ of the nominal value for D*^ (~ 3a) and within 
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±1.5 MeV/c^ for D*+ (~ Scr). 



The D and D* candidates are combined witli tiie K^Kg and K^K*^ can- 
didates to form B~ and S° candidates. Two kinematic variables are used to 
extract the B meson signal, the energy difference AE — — E^^^ and the 
beam-constrained mass Mbc = ^JiE^^^"^ — {p'^Y, where E'^ and are 
the center of mass (CM) energy and momentum of the B candidate and i?beam 
is the CM beam energy. The intervals Mbc > 5.2 GeV/c^ and \l^E\ < 0.2 GeV 
are selected. The MC estimate of the Mbc resolution is in the range (2.7- 
3.0)MeV/c^; the A£' resolution ranges between 12.5 and 14.5 MeV for final 
states including 7r°'s and between 8.5 and 10.5 MeV for other final states. The 
background from B D^*^D~ decay modes, with D~ K~K^^*\ is removed 
by requiring \M{D-) - M{K-K^^*'>)\ > 20MeV/c2 (~ Aa). Only one B me- 
son candidate per event is accepted. In cases of multiple entries, the chosen 
candidate is the one with the largest likelihood value, which is calculated using 
information about the differences between nominal and measured D, D* and 
7r° masses and the charged particle identification. 

The continuum background under the T(4S) signal is suppressed by topolog- 
ical cuts that were optimized using MC to model the signal and data in the 
B mass sideband (5.2 < Mbc < 5.26) GeV/c^ to model background. The 
ratio of the second to the zeroth Fox-Wolfram moments [9] is required to be 
less than 0.5. The angle 6'^^^ in the CM between the thrust axes of the parti- 
cles forming the B candidate and all other particles in the event must satisfy 
\cos{9^f^J.)\ < 0.85. The angle 9% in the CM between the beam direction and 
the B momentum direction should he in the range |cos(^^)| < 0.9. 



4 Systematics and Backgrounds 

The systematic uncertainties described below were estimated and added in 
quadrature to obtain final systematic errors separately for each decay channel. 
The contributions are listed in Table 1. 

A significant systematic uncertainty comes from the uncertainty of the charged 
track reconstruction efficiency, which has been evaluated from data. This un- 
certainty is estimated to be 2% per track for tracks with momentum larger 
than 200 MeV/c. A special study was performed to determine the momentum 
dependence of the reconstruction efficiency of the charged low momentum pion 
{P^^iTTsiow) < (220 - 250)MeV/c) from the D*+ meson decay. The distri- 
bution of the helicity angle, which is strongly correlated with the slow pion 
momentum, is compared in MC and data for B^ 0*^71" decays. Good 
agreement of the shapes of these distributions is obtained. A similar method 
is used to estimate the reconstruction efficiency of a 7r° meson produced in D*^ 
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decay. A test of the relative reconstruction efficiencies for low momentum tt"*" 
and 7r° was performed by comparing two inclusive D*^ decay modes. The ratio 
of the branching fractions Rd*+ = Br{D*^ — > D~^t:°)/ Br{D*+ — > D°7r+) was 
determined in generic BB MC and data. Reasonable agreement was obtained. 
A 10% systematic uncertainty is added to take into account the uncertainty 
in the reconstruction efficiency for both charged and neutral slow pions. 

The charged kaon and pion identification efficiency is checked by fitting the 
D*+ ^ D°n+, D° K-7r+ signal in MC and data. The fractions of signal 
events with and without kaon or pion identification are compared in the D*'^ 
center-of-mass momentum range from 1.0 to 2.5 GeV/c. This systematic error 
is found to be around 1% per particle. A similar procedure is applied to check 
the Kg selection efficiency, where the D^tt"*", —>■ /T^tt+tt" decay 

channel is used. The systematic uncertainty obtained is 3%. 

The non-resonant background under the K*^ signal is estimated from a fit 
using a relativistic Breit-Wigner function for the signal and a linear function 
to describe the non-resonant background. The mass and width of the K*^ are 
in good agreement with the nominal values. The non-resonant contribution 
under the K*'^ is found to be (2 ± 2) %. No correction is applied and an 
uncertainty of 4% was added for this background source. The efficiency of the 
D meson mass and vertex constrained kinematic fit is close to 100% with about 
a 1% uncertainty. The uncertainties in the D and D* meson decay branching 
fractions were also taken into account [5]. 

The probability of finding more than one B candidate per event is generally 
not large and is well reproduced by MC. The fraction of multiple entries varies 
from 2% for channels with only charged tracks in the final state, to 17% for 
channels with a D*^. The uncertainty in the signal yield due to this effect is 
estimated to be 4% for the decay modes with a D*'^ and 2% for other decay 
modes. 

The background under the B signal is separated into the combinatorial back- 
ground, which is removed by the fit procedure, and peaking backgrounds, 
where the content of the final particles is the same as in the studied channel, 
but a D meson was not formed. The systematic uncertainty of the combinato- 
rial background subtraction is estimated by varying the function that describes 
the background shape. This uncertainty is smaller than 5%. The width of the 
signal Gaussian is also varied inside the acceptable range and this systematic 
error is found to be in the range (3-5)%. The uncertainty in the estimate of the 
peaking background is found to be approximately 5% for the B~ D^K~K^ 
and B~ D^K~K*^ decay channels and significantly smaller for other chan- 
nels, if this background is evaluated from D meson sideband studies. 

The MC distributions of K~K*^ and K~Kg mass were adjusted to describe 
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the data. The shapes of these distributions are varied within their experi- 
mental uncertainties. For modes with Z)*° and D*^ mesons the polarization 
is varied from longitudinal to transverse. The uncertainty in the determina- 
tion of the number of BB pairs for the full data sample is 1%. Finally, the 
MC statistical uncertainties are also included in the systematic uncertainty. 
The overall systematic uncertainty ranges from 15 to 19% for modes with D 
mesons and from 19 to 22% for modes with D* mesons (Table 1). 

Table 1 

The systematic uncertainties in the branching fraction calculations. 



1. 


Reconstruction efficiency 






a) Charged track (per track) 


2% 




b) Slow 7r+ from D*+ 


10% 




c) Slow vr'J from D*^ 


10% 




d) K^, K*^, and mesons 


(1-3)% 


2. 


TT /K identiiication (per particle) 


1% 


3. 


Background subtraction 






a) Combinatorial background shape 


(1-5)% 




b) B signal width 


(3-5)% 




c) Combinatorial background under 


(1-5)% 




d) Non-resonant background under K*^ 


4% 


4. 


MC model 






a) D*° and D*+ polarization 


(4-8)% 




b) K^K^^*^ mass shape 


(2-3)% 




c) B candidate multiple entries 


(2-4)% 




d) MC statistical error 


(3-6)% 


5. 


D and D* branching fractions 


(2-7)% 


6. 


Number of B mesons 


1% 




Total: D modes 


~ (15 - 19)% 




Total: D* modes 


~ (19 - 22)% 



5 Results 

After applying all the selection requirements, the signal yields for the eight 
decay channels arc obtained from fits to the A£' distributions shown in Fig. 2. 
Because the /S.E and M^c parameters are almost independent, an additional 
selection Mt^c > 5.272 GeV/c^ was applied to suppress background. The B 
signal is described by a Gaussian function with a width fixed from MC. The 
peak position is fixed to zero if the signal statistical significance is less than 
4o". For decay modes with large numbers of events, the peak position is allowed 
to float. The background is described by a first-order polynomial, which fits 
events in the B meson mass sideband (5.2 < Mbc < 5.26)GeV/c^, shown 
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as the hatched histograms in Fig. 2. The backgrounds under the signals are 
well reproduced by the B mass sideband events. The fit interval —0.12 < 
AE < 0.2 GeV is chosen to avoid reflections from modes with an additional 
soft pion, in particular from D* mesons. The Mbc distribution for the range 
|A£^| < 25MeV was fitted by a Gaussian function with a fixed width to 
describe the signal and by the so-called ARGUS background function [10] with 
floating exponential and normalization parameters to describe the background. 
The peak position is fixed to the value 5.2795 GeV/c^ if the signal statistical 
significance is less than Aa. The shape of the background is chosen to be fiat 
if the population outside the signal region is small. The signal yields obtained 
from the AE and Mbc fits are listed in Table 2 and found to be in good 
agreement. 

Table 2 

The signal yields (obtained from the AE' and M^c fits), efficiencies (for Ktt and 
i^vrvrvr decay modes), branching fractions, upper Hmits and significances (the prob- 
abilities of a background fluctuation, measured in equivalent a units) for the studied 
decay modes. The branching fractions are calculated using the yield from the AE 
fit. 



Decay modes 


Yield 
AE 1 


Eff.(%) 

KTr/KTrTTTT 


Br. fractions, 
upper limits (10""^) 


Signif. 
a 




46.7 ± 8.2 

iG.-i ±7.4 


7.71 


7.5 ± 1.3 ± 1.1 


8.0 


BO ^ D+K-K*^ 


87.7 ± 11.4 

88.8 ± 10.2 


5.23 


8.8 ± 1.1 ± 1.5 


10.4 


B- D*^K-K*^ 


32.8 ± 7.2 
37.3 ± 6.9 


2.72 
0.99 


15.3 ± 3.1 ± 2.9 


6.7 


BO D*+K-K*^ 


37.5 ± 6.4 

38.6 ± 6.3 


3.28 
1.05 


12.9 ± 2.2 ± 2.5 


9.5 


B- D^K-K^ 


23.7 ± 5.9 
28.1 ± 5.8 


10.25 


5.5 ± 1.4 ± 0.8 


5.5 


BO D^K-K^ 


10.3 ± 5.0 
5.8 ± 4.5 


6.62 


1.6 ± 0.8 ± 0.3 

< 3.1 (90% CL) 


2.6 


B- D*°K-K° 


9.1 ± 3.9 
10.5 ± 3.5 


3.36 
1.26 


5.2 ± 2.7 ± 1.2 
< 10.6 (90% CL) 


2.5 


BO D*+K-K° 


5.4 ± 2.5 
5.6 ± 2.7 


4.46 
1.49 


2.0 ± 1.5 ± 0.4 
< 4.7 (90% CL) 


2.5 



Significant signals are observed in the channels with K K*^ subsystems and 
in the — > D^K~K^ decay mode (Table 2). The signals in the B^ 
D+K-K% B- D*^K-K% and B^ D*+K-K% decay modes are not 
statistically significant, and 90% confidence level (CL) upper limits for their 
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D^K k" 







\ d) 1 






jMjLh^ II ^ ii^i ini4 




0.1 0.2 -0.2 -0.1 

A E (GeV) 



Fig. 2. The AE distributions for eight B meson decay channels (see Table 1). The 
points with error bars are the data, the curves show the results of the fits described 
in the text. The hatched histogram shows the AE' distribution for events in the B 
mass sideband (5.2 < Mbc < 5.26)GeV/c2. 

branching fractions are also given in Table 2. We calculate the upper limits 
assuming a Gaussian distribution of the statistical error, and then inflate the 
limit by one unit of the systematic error. The branching fractions for the 
decays of the intermediate D, D*, and K*^ states were not included in 
the efficiencies and their values were taken from Ref. [5]. An equal production 
rate for the neutral and charged B mesons is assumed in the calculation of 
the branching fractions. 



For further analysis the events from the B meson signal region are selected 
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Fig. 3. a) The K~K*^ mass distribution after background subtraction and efficiency 
correction for the four B decay modes combined. The curve shows the fit to the 
relativistic Breit-Wigner function described in the text. The dashed histogram shows 
the distribution for the three-body phase space B decay, b) The shape of the ai 
resonance invariant mass for the a^(1260) p^ir" decay channel is shown by the 
sohd hue. The dashed hue (hatched area) shows the shape of the K~K*^ mass 
distribution for the a'^{\2QQ) —>■ K~K*^ decay channel. 

by the requirements \AE\ < 25MeV and M^c > 5.272 GeV/c^. The K'K*^ 
invariant mass distribution after background subtraction and efficiency cor- 
rection is shown in Fig. 3a for the combined signal region events for the four 
corresponding B decay modes. The efficiencies as a function of K^K*^ mass 
are flat within the errors, except for the bin where the D~ signal is removed. 
The background under the B signal is modeled using the Mbc sideband. 

Simulation of the three-body decay B D^*^ K~ K*^ according to phase space 
results in the K^K*^ mass distribution shown by the dashed histogram in 
Fig. 3a. It is clear that the observed shape and an enhancement in the low 
K*^ mass region rule out any significant phase space contribution. A MC 
simulation assuming an intermediate K~K*'^ resonance shows that the ob- 
served low mass enhancement can be explained by the quasi-two-body decays 
B 0^(1260), with 0^(1260) K" K*^ . Figure 3b shows the invariant 

mass of the products of the af (1260) meson decay via its dominant decay 
mode p^TX^ (the solid line) and via the decay mode K^K*^ (the hatched 
area). The shapes of the distributions were described by relativistic Breit- 
Wigner functions where both decay channels contribute to the total width 
simultaneously [11]. 

The K~K*^ mass distribution was parameterized by the sum of the Breit- 
Wigner function for a resonance and a phase space function Fps{Mkk*) with 
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a fixed shape for the non-resonant three-body phase space B decay, 



dMKK* {Ml - Ml^.Y + Ml TUMkk*) 



+ fps ■ Fps{Mkk*) , (1) 



where the mass dependence of the partial TpartiMKK*) and total TtotiMxR*) 
decay widths is approximated assuming a two-body phase space oi — > KK* 
decay: 

^part{MKK*) = ^KK* ' {QKR* / Qo) , Ttot{MKK*) = ^ai ' (?p7r/?o)- (2) 



Here /^^ and fps are normalization parameters, and M^^ and are the mass 
and the width of the ai resonance. The function Fps{Mkk*) describing the 
three-body phase space distribution obtained by simulation was parameterized 
by a product of polynomial and threshold functions. The QpT^, Qkk* and 
parameters are three- vector momenta of each daughter particle in the ai rest 
frame. The constant go is calculated assuming a ai — p7r decay with the mass 
of the fli resonance fixed to its nominal value. The total width TtotiMxR*) is 
approximated by the dominant partial width of the ai — > pvr decay mode. 

If the ai mass is fixed to its nominal value Ma^ = 1230MeV/c^ [5] and the 
width is fixed to any value in the range 300 < < 600 MeV/c^, the relative 
contribution of the ai extracted from the fit is almost 100%. Introducing an ad- 
ditional interference term in the formula (1) leads to a small correction, which 
does not change the fit results. If the phase space term is omitted, a fit with 
the resonance width fioating gives Ta^ = (460 ± 85) MeV/c^ (Fig- 3a), which 
agrees well with the nominal range for the ai width [5]. The quality of the 
fit is reasonable: /n.d.f. = 29.5/18. Unfortunately, contributions from non- 
resonant three-body B decays that do not have mass distribution with a phase 
space shape and/or higher resonances cannot be reliably separated due to the 
low statistics. The branching fractions for B D^*)a7(1260) have been mea- 
sured by CLEO in the Oi fnr mode [3] . Comparing the CLEO measurement 
with our results and assuming purely resonant production of the K~K*^ sys- 
tem, the relative fraction R = 5r(ar(1260) ^ R- K*^) / Br{a^{l2m) all) 
is found to be within the range (8-15)% depending on the B decay mode 
studied. This is larger than the value (3.3 ± 0.5)% that CLEO obtained by 
studying the three pion mass shape in r decays [4]. 

To study the quantum numbers of the intermediate K'K*^ state, the angular 
distributions of the final state particles were examined. First, the helicity angle 

9kk was defined as the angle between the momentum of the K~K*^ system in 
the B meson rest frame and the momentum of the K*^ in the K~ K*^ system 
rest frame. The helicity angle Ok*-, which was defined as the angle between 
the momentum of K*^ in the K~K*^ system rest frame and the momentum 
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COS 6 cos 



Fig. 4. The cos Ork distributions: a) for the sum of the B D^K K*^ 

and ^0 ^ D+K-K*^ modes and b) for the sum of the B' D*^K-K*^ 
and B^ — > D*'^K~K*^ modes. The cos 9k* distributions: c) for the sum of the 
B- D^K-K*^ and 5° ^ D+K'K*^ modes and d) for the sum of the 
B- D*^K-K*^ and 5° -> D*+K-K*^ modes. The curves are fits to a con- 
stant (dashed), sin^^ (dotted) and cos^^ (dashed-dotted) functions. 

of the daughter K'^ meson in the K*'^ rest frame, was also studied. 

Table 3 gives the expected distributions for these two helicity angles for the 
possible quantum numbers of the K~K*^ system in i? — > DK~K*^ and 
B — > D*K~K*° decays. One can see that for = 0~ and 1+ the angular 
distributions for the modes with D and D* are the same. However, for modes 
with a D* and = 1~ the situation becomes more complicated and the an- 
gular distributions depend on the fraction of the longitudinal polarization in 
the B decay. 

Figure 4 shows the cos 9kk and cos 9k* distributions for the first four B decay 
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Tabic 3 

The angular distributions of the K^K*^ system with = 0^, 1^, 1^. The values 
of /n.d.f. were obtained from fits to the experimental angular distributions (see 
Fig. 4). 





D meson modes 


D* meson modes 


Sum 
X^/n.d.f. 


Okk 


Ok* 


Okk 


9k* 


0- 


const 


cos^ 6 k* 


const 


cos^ 6 k* 


71.7 / 16 


1- 




SIIL^Ok* 






37.3 / 8 


1+ 


const 


const 


const 


const 


20.7 / 16 



mode combinations after background subtraction and efficiency correction. 
The angular resolution in these measurements is significantly smaller than the 
bin size. No significant deviation from the fiat dependence is observed for all 
studied angular distributions, except for the bin at high cos 9kk in Fig. 4b. 
The results of the fits are also shown in Fig. 4. The ■x^/n.d.f. values given 
in the last column of Tabic 3 were obtained by summing the values for the 
four individual fits. The fit results are consistent with a pure J^=l+ state for 
the K'K*'^ system and rule out pure J^=0~ or J^=l~ states. It has to be 
stressed that only decay modes with D mesons are used to test the angular 
distributions for J^=l~ case. The fit results support the hypothesis that the 
intermediate aj^(1260) resonant state decays to the K~K*^ final state. 



In the four decay modes to K~Kg final states, a clear signal is only ob- 
served in the B~ D^K~K^ decay mode (Fig. 2e), the enhancements in 
the other three decay modes (Figs. 2f,g,h) each have less than 3o" significance. 
The K~Kg mass spectrum for the B~ — > D^K~ decay mode after back- 
ground subtraction and efficiency correction is shown in Fig. 5a, where the 
— > ir~7r+ decay channel is used. The mass distribution peaks near the 

> D^K~K^ signal events in the 



K 



'K^ mass threshold; the fraction of B 



mass range M{K K^) < 1.3GeV/c^ is 55%. The distribution for the three- 
body phase space B decay is shown by the dashed histogram. 



The efficiency corrected and background subtracted heficity angle distribution 

for the K~Kg system is shown in Fig. 5b. The angle 6kk is defined similarly 
to that for the K~K*^ system, where the is substituted for the K*^. This 
distribution is fitted to the form. 



dN 



dcos Okk 



A R 



L COS Ukk 



[I - Rl) sin^OKK 



(3) 



where Rl = F^/F is the fraction of longitudinal polarization and A is a 
free normalization parameter. A fit to the functional form (3) gives the value 
Rl = 0.97 ± 0.08, which is typical of a P-wave decay of a =1~ system to 
two pseudoscalar mesons. 
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M(K" kJ) (Ge\/c) cos Q 



Fig. 5. The fully corrected a) K~Kg invariant mass and b) K~Kg helicity angle dis- 
tributions. The data are from the B~ D^K~Kg decay mode with K~tt~^ . 
The dashed histogram shows the expected distribution for three-body phase-space 
B decays. The curve shows the result of the fit to the function described in the text. 



The observed mass distribution is expected to include interfering contributions 
from the p(770), p(1450) and p(1700) resonances and follow that for the KK 
system produced in the vector isovector state in other processes like e+e^ anni- 
hilation into two kaons or t~ lepton decays to K~K^. Although rather precise 
measurements of the cross sections exist for the processes e^e~ K^K~ [12] 
and e+e~ — ^ K^Kg [13], the final states observed in e+e^ annihilation are a 
mixture of isovector and isoscalar contributions and their separation requires 
a complicated mo del- dependent analysis, r lepton decays are free of this dis- 
advantage since in this case the produced K'K° system always possesses 
the necessary quantum numbers. CLEO [14] studied this mass distribution 
in T decays and concluded that the observed K^K^ mass distribution could 
qualitatively be described by a p-like intermediate mechanism producing an 
enhancement close to the K^Kg threshold. 

The measured branching fractions for the decay modes with the K~Kg system 
are of the order of (2—4) % relative to those of the respective B D^*^ p~ (770) 
decay modes with p^(770) 7r~7r°. These values are at least 5 times larger 
than one would naively estimate, using r decay branching fractions corrected 
by phase space factors [15], even taking into account positive interference of 
the p(770) and p(1450) resonances. 
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6 Conclusions 



In summary, eight B — > decay modes have been examined for 

the first time and significant signals have been found for five of them. The 
branching fractions and upper Umits are obtained. Significant signals are ob- 
served for the four modes with a K^K*^ system. The K*^ mass spectrum 
is peaked near threshold and smoothly decreases with increasing K~K*^ mass. 
The angular dependences indicate that the K''K*^ system has J^=l+. The ob- 
served behavior can be interpreted as production of an intermediate aj"(1260) 
resonance that then decays to K~K*^. 

A clear signal is also found in the B~ — > D^K'Kg decay mode. The K~Kg 
mass spectrum for this signal is also peaked near the K~Kg mass threshold. 
The angular dependence shows a cos^Okk behavior indicating that the K~Kg 
system has J^=l~. Some enhancements in the signal regions are also seen in 
the B- D*^K-Kl, B^ D+K-R^g and B^ D^+K-R^ decay modes, 
but with limited statistical significance. A non-resonant three-body contribu- 
tion, which is not described by phase space, cannot easily be ruled out for any 
of these decay modes. 
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